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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGATION  OF 
THE  EFFECT  OF  DYES  ON  SOLAR  DISTILLATION 


By 


Cochairman:  C.  K.  Hsich 

Major  Department:  Mechanical  Engineoring 

An  analytical  and  experimental  study  of  the  effect  of  dyes  on 
solar  distillation  was  conducted  in  this  study.  The  analytical  model 
developed  treated  the  transient  heat  transfer  inside  the  dye-water 
system  as  one  dimensional.  The  bulk  fluid  was  discretlsed  into  layers 
with  conduction,  convection  and  radiation  interactions  occurring 
between  them.  Coupled  with  this  was  the  energy  exchange  from  the  sur- 
face of  water  to  the  glass  cover  by  convection,  radiation  and  evapora- 

insulated  and  a black  absorber.  A finite  difforencc  technique  was  used 
in  solving  the  nonlinear  partial  differential  equations.  In  the 
computer  program  the  measured  solar  radiation,  ambient  temperature  and 

te  history  of  the  water  layers  and  glass  cover  and  the 


distiller 


solar  stills  (4  ft.  x 4 ft.  x 1 ft.).  The  control  still  contained 
water  only  while  the  other  had  dye-water  solution.  The  dyes  used  were 

trations.  The  Green  dye  was  made  by  mixing  a 33  percent  by-weight 
mixture  of  Neptune  Blue,  Tartraxine  and  Red  Carmolsinc  (all  GAP  com- 

by  a specific  dye  over  that  from  the  control  unit  was  developed.  The 

for  clear  days.  However,  no  difference  in  the  still  productivities, 
between  the  control  and  with  dye,  was  noticed  on  a completely  cloudy 
day.  An  empirical  relationship  was  developed  for  evaporation  heat 
transfer  coefficient  and  the  distillate  output  as  a function  of  brine 

suitable  from  two  points  of  view:  Increased  evaporation  and  lack  of 

duction  of  the  absorption  coefficient  by  as  much  as  95  percent  at  0.5  un 
wavelength  after  a month  of  exposure.  The  degradation  of  Green  dye  was 
not  as  severe  as  the  Red  dye,  and  tills  was  primarily  because  of  the  Red 


excellent. 


s of  ambient  temperature,  wind  velocity  and  dye  concentration  on 

case  in  output  with  increase  in  ambient  temperature  and  an 
se  in  still  productivity  with  wind  speed. 

output  resulted  with  increasing  dye 

lose  tested)  with  Black  having 


different  for  different  dyes  (i 


distillate  with  concentration  times  concentration  and  the  distillate 

Gainesville,  Florida,  the  optimum  concentrations  for  Black,  Green  and 
Red  dyes  were  218  ppm,  377  ppm  and  d08  ppm, 
cions  for  further  research  were  given. 


CHAPTER 


INTRODUCTION 


Abundant  it  is,  but  not  infinite.  Of  the  453 
; 10“  gallons)  of  sea  water  that  evaporates  every 
y 38  units  are  available  for  possible  uses  by 
se  38  units  only  70  percent  can  be  utilised  11] 

enough  supply  of  fresh  water  and  the  present 
that  the  hydrological  cycle  can  support  about  25 


le  available  annual  water  supply  will  probably  be  insufficient  f< 


.n  conjunction  with  the 
rater  shortage.  Already 


by  transporting  fi 


drinking.  Even  in  areas  well  endowed 


with  ground  water  supply,  because  of  increased  Industrialization  and 
population  growth,  the  cate  of  depletion  of  these  supplies  is  reaching 
alarming  proportions.  With  the  cost  of  water  transportation  already 


.s  problem  of  water  shortage 


high,  some  other  way  of  getting  fi 

about  7 .5  x 10  gallons  of  water  per  day  (gpd)  in  about  a dozen  coun- 
tries [3].  Among  the  processes  used  about  85  percent  are  Multistage 

technology  in  this  area.  The  energy  used  in  these  plants  is  either 


However,  with  Increasing  fossil  fuel  shortages  and  price  r. 
is  but  natural  to  look  at  some  other  methods  ol 


heating 


facility. 
Arabia)  bas  bean 
, nevertheless,  i 


plant  (capacity  of  5,000  gpd  in  Saudi 
using  the  second  approach  [51 . It  is. 

average  of  about  500  gpd  of  water  is  used  for  drinking  and  cooking 
purposes.  In  these  communities  setting  up  a large  solar  desalination 

transportation  costs.  Several  research  workers  have  therefore  built, 
and  experimented  on,  medium  sire  solar  distillation  units  (output  of 

A solar  distillation  unit  (solar  still)  is  a very  simple  device  in 

present  chapter  it  is  sufficient  to  describe  briefly  the  working  of 
such  a still  which  is  schematically  shown  in  Figure  1.  Saline  or 

having  depth  ranging  from  1/2  inch  to  approximately  1 to  3 ft. 

shallow  basin  (depth  - 2 inches)  or  deep  basin.  The  bottom  of  the 


nc,  continuously  or  periodically.  Above  tl 


research  workers  all  over  Che  world  bul 


i,  however,  choc  mosc  of  chese  variacions  in  design 
drawbacks  which  have  made  them  cechnically  and 

discillacion  unics.  Thus  chere  is  a need 

< an  attempt  co  cesc  Che  feasibility  of  a novel  concept 

tee,  thereby  increasing  Che  productivity  of  Che  still. 


This  survey  covered  both  the  large  installations  (output  of  < 


thereby  reducing  the  productivity  of  the  still, 
the  high  initial  cost,  which  made  them  unattractive  for  other  communi- 


e for  large  or  even  small  t 
for  a better  distillation  unit.  But 

ductivity  of  solar  stills.  In  addit 

the  performance  of  a still  are: 


experimentally  investigated . He 


study  will  alt 


in  Figure  2.  The  data  in  this  plot  are  iron 


most  of  the  distillation  output  of  shallow 


of  solar  radiation  is  absorbed  in  the  water  and  because  of  the  thermal 


cycle.  In  addition. 


Scientists  at  the  University  of  California  have  studied  the 
cascade  type  still  [10].  The  depth  of  the  brine  was  about  1/4  inch. 
High  output  rates  were  obtained  (0.13  gal/ft2-day  at  about  2000  Btu/ft2 


still  was  beset  with  problems  like  dry 
of  black  epoxy  paint,  deposition  of  salt 


l,  blistering 


Distillate  (pnllons/ft^-day) 


is  required  expensive  maintenance,  t 


s til)  hi 


tills,  A porous  black  material  s 
ne  water.  Since  the  layer  ot  tin 
s rapid  evaporation  and  increase 
0.16  gal/ft  -day  at  about  2000  Bl 


if  productivity  results  (out- 
pores,  dry  spots  developing 


resulting  in  deterioration  of  the  wick  material  and  precise  contrc 

attractive  for  commercial  production.  Recently  ESB  Corporation  |J 
patented  a new  type  of  wick  material  which  is  claimed  to  have  se] 


cleaning  properties 
detailed  data 

Several  investigators  In 
Recently  Szulmayer  [16]  has  reported 


lot  clogged  b; 
difficult  to 


similar  shallow  basin  stil 


tr  Material  ai 


lover  material  has  two  functions.  Besides  allowing  solar 


>f  condensing  surface  ai 


radiation 


n (15)  tl 


i of  the  glass,  more  than  offsets  the  increa: 

Normally,  two  types  of  cover  materials  are 
advantages  of  glass  are: 


high  wettability  of  glass  and  relatively  high  stability  of  properties 

requires  an  increase  in  its  thickness  (normally  0.12S  inch  thick  glass 
is  used  for  cover) . Plastic  films  normally  used  for  solar  stills  have 
high  solar  transmission  (an  average  of  90  percent  for  0.004  inch  thick 

pared  to  about  85  percent  transmittance  for  0. 125  inch  glass) . At  the 
same  time  these  films  have  high  Infrared  transmittance  (80  percent  in 
Che  wavelength  region  of  4 pm  to  7 pm)  and  high  susceptibility  to  ultra- 

of  plastic  cover  unattractive  in  solar  stills.  Thus,  till  the  advent 
of  better  plastic  material  which  will  overcome  the  difficulties  listed 


the  whole  day,  and  c)  tl 

both  experimentally  (16] 

solar  radiation  entering  the  still  is  small,  b 

the  major  factor  governing  the  shape  of  the  co1 
of  enough  condensing  area.  It  is  interesting 
investigators  appear  Co  have  considered  this  f; 

and  the  water  surface  which  ultimately  results 


it  intercepts  maximum  S( 
available  then  the  effect  of  spacing  bel 


Ambient  Temperature  ai 


id  Velocity 


These  two  variables  tu 


the  productivity. 


coefficient  thereby  increasing  the  heat  loss  from  the  cover  and  the 
productivity.  However,  Lof  cc  al.  (18),  Baum  [19]  and  Bloomer  ct  al. 
[20)  have  predicted  that  increase  of  wind  velocity  decreases  the  output . 

is  not  borne  out  by  experiments  (21),  calculations  [22J  and  by  computer 


ie  condensing  surface  1; 


i function 


condensing  the 


productivity  of 
investigators  have  tried  t( 

investigated  a multiple  effect  solar  still  where  the  heat  of  condens. 

in  the  second  effect  and  so  on.  The  condenser  of  the  last  effect  wai 
cooled  by  water.  Salam  and  Daniels  [25)  used  two  inflated  plastic  tt 


external  condenser  where  t 

between  20  and  40  percent 
efficiencies  this  system  : 

Although  production  i 

equipment  for  air  blowing 


2 gal/ft  -day  on 
the  additional  cc 


e also  reached  by  Dunklc  [24J, 

vapor  tightness;  no  or  minimal  seepage  of  brine  from  the  basin;  good 
frequency  of  flushing  and  the  orientation  of  the  still  ( 26 ] . 


ic  above  features.  By  mixing  water  soluble  dyes  ii 


surface  layer 


advantage  of  cooler  condenser  surface.  The  problem  of  salt  buildup 
bottom  of  the  basin  and  does  not  interfere  in  the  solar  radiation 


: of  2-Napthol  Green  dye  on  thl 
ir  different  c« 


cent  for  brine  layers  of  20  cm  and  50  cm  deep) . They  did  report 

19  percent  over  that  of  uncolored  brine  was  achieved.  However,  no 
empt  was  made  to  study  the  effect  of  other  dyes,  the  effect  of  sun- 
ht  on  the  degradation  of  the  dye  or  to  correlate  their  experiments 

ame  on  accepted  practice.  However,  very  little  data  are  available 
cess.  The  use  of  red  bacteria  Halobaoterium  or  blue  green  algae 


Spirulina  has  been  tried  to  enhance  the  evaporation  of  brine  I 29 ] . The 
advantages  of  such  a system  are  obvious.  There  is  no  solar  degradation 
of  the  bacteria  and  at  the  same  time  the  culture  can  be  obtained  free 

costly  dyes  altogether,  which  are  not  anyway  recoverable.  However, 

tion  and  temperature  and  have  not  been  therefore  commercially 
successful  [30). 

Keyes  et  al.  [31]  have  extended  Block's  work  to  include  other  dyes 
and  tested  their  degradation  by  sunlight.  Among  the  six  dyes  tested 


‘ salt  solutions  and 
conflicting  reports 
in  solar  evaporation. 


The  present  investigation  is  an  attempt  to  achieve  the  following 
objectives: 

These  parameters  are  ambient  temperature,  wind  velocity,  type 
prediction  of  temperature-time  history  of  the  dye-water 
ditions  and  input  solar  radiation.  The  model  is  developed  in 

productivity  of  deep  basin  solar  still.  The  productivity 
from  this  still  is  compared  with  that  from  an  identical 

used  in  this  study.  The  dyes  used  are  Napthylamine  (black), 
e (Dark  Green) . 


CHAPTER  III 
THEORETICAL  ANALYSIS 


In  this  chapter  a theoretical  analysis  for 


Physical  Model 

analyxe  Che  energy  transfer  process  in  the  basin  the  dye-water  system 
has  been  divided  into  discrete  layers  of  thickness  As.  An  energy 

assumptions: 

e large  compared  to  those  in  the  x direction 


s,  the  changes  in  physical  properties  ol 


•ffl— 


solution  like  viscosity,  density,  specific  heat  and  thermal 
conductivity  are  negligible. 

The  spectral  emission  and  scattering  coefficient  of  the  dye- 


with  condensate  f: 
averaging  out  the 
respect  to  angle  of  incidence. 


le  above  assumptions,  the  energy  balance 


similarly 


* "3 

— —,•§=  £■=?—] 


*N' 

1 

™--^Att=Sr]  (W 

M •[--]■  §=.-®=^ 

(7b) 


surface  is  considered  as  having  wavelengths  in  the  range  of  0.3  inn  ■ 
2.0  pm.  The  choice  of  this  range  has  been  dictated  by  the  fact  tlia 


he  glass  surface.  Furthermore,  to  facilitate  computation,  ( 
.r  spectrum  132]  has  been  smoothed  by  fitting  a polynomial  i: 


_ j m f 23791  - 713.7;  0.2  urn  < 
“ 1 1261.2  exp(-1.951);  0.5 


modified  for  different  times  * 


s given  by  q(t>  then  q^(t)  h 


eft) 


Absorption  of  Radiation 

The  incident  angle  of  solar  radiation  on  the  dye-water  surface 
incidence  changes  the  absorption  thickness  for  any  layer  of  thickness 


index  of  refraction  for  water-dye  mixture. 


compared  to  the  air  inside  the 


From  Fresnel  equation  the  reflectance  can  be  calculated  by  knowing 


x sln^Oj  - er) 


(13a) 


dye-* 


where  q^Ct)  is  obtained  from  equation  (9),  and  r is  given  by  equations 
(12)  and  (13)  above. 


4x*  = 4x/cos0c 


(14) 


q"  jtOlexpf-ajAx*)  )d^ 


(15) 


is  the  spectral  absorption  coefficient  of  the  medium.  The  absorption 
coefficients  of  different  dyes  have  been  obtained  experimentally  and 
will  be  discussed  in  greater  detail  in  Chapter  V. 


28 


where  t = thickness  of  layer 

The  data  for  equivalent  conductivity  for  plane  layers  from  different 
investigations  has  been  plotted  134]  and  fall  on  a straight  line,  given 


l°SlO 


where  k ■ conductivity  of  the  water 

KS£  » Raleigh  number  based  on  thickness  t of  the  layer 
„ Bg  (TH  - T.K3 


diffusivity 


Convection  Li 


interferometric  studies  on  a laboratory  still  and  found  that  most  of 

enclosure.  The  bulk  of  air  at  the  center  of  the  enclosure  remains 
essentially  stationary.  Thus,  the  evaporation  rate  is  determined  by 
Che  analogy  between  heat  and  mass  transfer  for  free  convection. 

The  heat  transfer  by  convection  from  Che  watur  surface  to  the 


given  by: 


simplicity  an  average  temperature  between  the  top  and  bottom  port i< 
The  heat  transfer  coefficient  h can  be  calculated  from  the 


the  water  surface  and  the  glass  cover)  for  most  stills.  Experimentally 

The  bouyancy  factor  (dp/p)  in  the  Grashoff  number  of  equation  (20) 
vecting  fluid.  However,  this  has  to  be  modified  for  the  solar  still 


since  in  this  case  the  humidity  of  air  changes  along  the  fiow  path. 
This  modification  is  achieved  by  assuming  that  air-water  vapor  mixture 
is  an  ideal  gas  [36] . 


Therefore,  equation  (21)  is  modified  ti 


Thus,  the  density  of  water-air  mixture  near  the  glass  cover  and  at  the 
water  surface  can  be  respectively  written  as: 


(23) 


32 


h - C t ^ (S(Ap/p)  U CP)  1/3 


' mass  transfer  coefficient,  mass/arca  - til 
driving  potential  for  the  mass  transfer. 


dimensionless 


m j ■ __i  ■ mass  fraction  of  species  i;  (where  subscript  s implies 


Prom  analogy  bt 


transfer,  an  equation  similar 


Cr  • Grashoff  number  for  mass  transfer 
Num  - mass  transfer  Kusselt  number  ■ m L 


■ binary  diffusion  coefficient  o. 
* Schmidt  number  = Vci 


dif fusivity 


v * kinematic  viscosity  of  the  mixture 


The  above  analogy  is  based  on  the  following  assumptions: 

(iii)  Negligible  viscous  dissipation 
(iv)  Negligible  emission  or  absorption  of  radiant  energy 
In  a solar  still  there  is  combined  heat  and  mass  transfer  from  the 

functions  of  Pr,  Sc,  and  the  ratio  of  bouyancy  induced  by  the 
temperature  difference  to  that  Induced  by  the  mass  concentration  dif- 


air.  Ulus  the  problem  is  simplified  in  the 


"ig*  hfg 


a transfer  coefficient 


(m,  -m.  )»  mass  transfer  driving  potential 

thus,  from  equations  (24)  and  (28)  the  mass  fraction  of  species  is: 


Similarly, 


The  difference  equations  have  been  written  for  different  sections 
and  have  been  thus  identified. 


(expt-OjAx*)  ] di) 


lte  • apparent  conductivity  given  by  equation  (16) 
equal  to  k.  The  boundary  conditions  can  be  written  as: 


Top  Laver 

Combining  equations  (3)  and  (4)  give  the  temperature  of  the  top 


emissivity  shape 


A ' <,!«*»  («*<-<■***>  ">  ■ 


outside  convection  heat  transfer  coefficient 

specific  heat  of  glass  and  aluminum  cover, 
respectively 

thickness  of  glass  and  aluminum  cover,  respectively 


■ other  symbols  a: 

Experimentally  it  was  : 


, - 15"F  [37 J 

ik  place  on  the  back  plate.  , 

temperature  was  equal  to  that  of  the  glass  cover.  Moreover,  i 
of  magnitude  analysis  of  equation  (37)  shows  that  the  storage 


38 


the  aluminum  is  negligible.  Thus,  all  the  temperatures  T . in  the 
above  equation  can  be  considered  as  T . This  yields 

Purine  Daytime 

OFeFA  It*  - T*)  + hmA  (mls  - mig)  hfg  + h.A(T*  - t')  + og  Ag  q(t) 


During  Nighttime 


(38) 


Method  of  Solution 


cop  and  bottom  surfaces.  Therefore,  a computer  program  was  written  to 

the  glass  cover,  simultaneously.  Explicit  method  of  solution  was  used 
in  solving  them  and  therefore  the  time  interval  At  was  carefully  chosen 
so  as  to  ovoid  unstable  solutions.  Appendix  IX  shows  the  stability 


present  analysis. 


getting  the  temperature-time  history  of  various  layers, 
initialized  at  sunrise,  using  their  experimental  values. 

ambient  temperature  and  wind  velocity, 
of  equation  (35)  is  by  regulafalsi 


3,  then  the  apparent  conductivity  k i 


calculated.  The 


(b)  All  the  layers  above  layer  n will  have  convection 
transport  with  AT  = Tg  ^ - Tg  and  subsequently  k is 
calculated  using  equation  (17)  and  depth  l . 

(c)  The  layers  below  the  nth  layer  have  heat  transfer  between 

Steps  2 to  4 are  then  repeated  with  an  increment  of  time  At, 
till  a 24  hour  cycle  is  completed. 


generate  temperature-time  history  for  different  layers, 
results  of  the  analytical  model  are  presented  in  Chapter  V. 


Temperature  of  layers  (arbitrary 


Tipure  6.  Temperature  profile  of  the  layers. 


CHAPTER  IV 

EXPERIMENTAL  SETUP  AND  PROCEDURES 


Distillation  Units 


troughs.  This  trough  also  provided  support  for  the  glass  cover.  A 


enhance  Che  distillate  fli 
using  Argon-arc  welding. 


i che  incerior  of  Che  basin.  The  glass 


a)  chere  should  be  sufficicnc  condensing  area  and  b) 


h facing  glass  C( 


s double  screngch  (3 


wich  hypotenuse  as  5 ft. 

solar  radiacion  as  well  as  condenser  surface. 

by  silicone  rubber  sealanc  Co  make  the  cover  vapor  tight.  The  vapoi 
tightness  of  Che  cover  was  Chen  checked  by  passing  water  from  a rubt 


ip.  This  compound  remains  pliable  over  cl 


n schematically  in  Figure  1- 


couples  were  attached  on  outside  of  the  Plexiglas  basin  by  a paste  of 
Plexiglas  chips  and  Ethylene  dichloride,  while  on  the  glass  cover  they 
were  held  down  firmly  by  transparent  scotch  tape.  This  tape  was  over 
an  aluminum  covered  styrofoam  piece  pressing  down  on  the  thermocouple 
as  shown  schematically  in  Figure  13.  Besides  tl 
couples  were  attached  on  Che  outside  of  the  insulation  ai 


f dye-water  in  the  stil 

in  Figure  14.  The  thermocouple  junctions  were  coated 
coating  of  shellac  to  make  them  corrosion  resistant. 
These  thermocouples  were  calibrated,  together  wi' 


There  were  about  30  thermocouples  attached  at  different  locations 
both  the  stills  and  this  necessitated  the  use  of  multiple  point 


recorders.  Two  such  recorders,  one  made  by  Honeywell  and  the  other  by 
Leeds  & Northrup,  were  used.  The  Honeywell-Brown  electronix  recorder 
had  20  points  with  accuracy  of  - 0.5°!'  and  range  of  0-150"F.  The  Leeds 

and  range  of  0-200°F.  These  recorders  were  connected  to  the  mechanical 
timer  so  that  the  temperature  readings  could  be  taken  after  every  1 

The  solar  radiation  data  was  measured  by  Eppley  Phyrheliometer, 
Model  Ho.  10,  and  recorded  on  a single  point  Leeds  & Northrup  speedomax 

temperature  readings  were  obtained  from  the  hygrothermograph  at  the 


,F  Corporation.  They  were: 

3.  Mixture  of  33  percent  by  weight  of  each  of  the  following 
Green  dye. 

The  absorption  spectrum  from  wavelengths  of  0.3  pm  to  0.75  pm  was 


(Incident  radiation/Transmitted  radiation}]  directly.  However,  the 
Perkin-Elmer  Monochromator  is  not  a differential  device  and  thus 

absorption  coefficient  o:,  was  then  calculated  from  these  readings. 


Figure 


tank.  Common  salt  (composition  unknown),  3.5  percent  by  weight,  was 
added  to  the  two  stills.  In  one  of  the  stills,  dye  of  known  concentra- 
tion was  added.  The  covers  were  then  put  over  the  stills  and  allowed 
to  stand  for  a day  or  two  to  attain  a steady  periodic  operation.  After 
this  period  of  time  the  distillate  was  collected  every  24  hours  and  the 
weight  of  the  distillate  was  measured  by  a cantilever  type  balance  with 
an  accuracy  of  i 0.1  lbs.  The  temperature  measurements  of  various 
locations  were  recorded  on  the  strip  chart  at  every  1 hour  interval. 

Three  dyes  were  tested  wich  concentrations  of  50  ppm  and  100  ppm 
for  each  except  for  the  Black  dye  where  the  testing  was  done  for  50  ppm 
and  172.5  ppm  concentrations.  The  duration  of  testing  for  each  dye  is 

On  a number  of  occasions  the  distillate  was  collected  every  hour 
so  as  to  study  its  hourly  variation.  This  also  helped  in  calculating 
Che  constant  C in  equation  (20),  as  will  be  shown  in  Chapter  V. 

The  absorption  spectrum  of  the  dyes  before  and  after  the  duration 

to  study  carry  over  of  the  dyes  by  water  vapor. 

The  results  of  the  experiment  are  presented  in  Chapter  V. 


The  maximum  concentration  of  Black  dye  was  supposed  to  be  150 
ppm  but  the  resulting  mixture  inadvertently  resulted  in  172.5  ppm 


Caraoislne  (Red) 


11/10/77  - 12/8/77 
12/10/77  - 02/14/78 
02/16/78  - 03/28-78 
03/29/78  - 04/24/78 
04/26/78  - 08/13/78 
08/15/78  - 09/21/78 


CHAPTER 


RESULTS  AND  DISCUSSIONS 


Dve-Absorp cion  Spectrum  Results 

The  absorption  spectrum  of  different:  dyes  is  shown  in  Figures  21 
through  25  where  the  absorption  coefficient  a,  is  plotted  against 
wavelength  X.  The  absorption  coefficient  is  a function  of  concentra 


present  analysis  the  temperature  and  pressure  effects 

(Red)  dye*  taken  both  before  and  after  the  experiment 
the  dip  in  the  absorption  curve  between  0.6  end  0.7  pm 
part  of  the  spectrum.  Similarly  in  Figures  22  and  23 


strictly  speaking,  is  a dark  blue  dye  ani 

Blue,  Tort ranine  (Yellow)  and  Carraoisine 
absorption  curve  in  the  visible  region  ii 
region  of  the  spectrum  (see  Figures  24  ai 


experimental  data.  ' A 


I 


the  dyes  tested  in  the  present  study  have  a,  very  nearly  the  same  as 

wavelength  region  most  of  the  energy  of  the  incoming  photon  goes  to  the 

water  complex.  Since  Che  amount  of  Che  dye  is  small  ('100  ppm),  the 

However,  in  the  wavelength  region  of  0.7  pm  to  0.9  pm  there  is  still 
enough  energy  (1.55  eV)  in  the  photon  to  produce  higher  by  changing 

dyes  tested  in  this  study  did  not  show  this  change,  and  thus  there  is  a 
need  to  develop  dyes  which  will  have  high  a.  in  this  region  since  it 

their  absorption  spectrum  before  and  after  the  experiment.  Figure  20 

95  percent  at  a wavelength  of  0.5  Mm.  This  also  resulted  in  the  change 

The  reason  for  this  degradation  appears  to  be  photolysis  caused  by 
ultraviolet  radiation  in  the  solar  spectrum  since  a part  of  this 
spectrum  (ultraviolet  region  from  0.3  pm  to  0.4  pm)  is  transmitted  by 


■Q- 


The  bond  strength  of  N=N  is  3.45  eV  [40]  while  the  energy  of  the 
photon  in  0.3  Mm  wavelength  is  4.13  eV  and  thus  the  degradation  of  dye 
is  inevitable.  However,  no  one  single  mechanism  explains  the  fading  of 
dyes  and  involves  three  photochemical  reactions — oxidation,  reduction 

photolysis  and  oxidation  reactions  are  possible  for  its  fading.  One 
possible  indication  is  the  change  in  pH  of  this  dye.  in  100  ppm  concen- 
tration solution,  from  5.2  (before  experiments)  to  6.7  (after  experi- 


the  Carmoisine  dye.  The  absorption  coefficient  of  this  dye  drops  by 
about  85  percent  at  0.5  Mm.  This  resulted  in  the  change  of  the  solu- 


tion color  from  dark  green  to  transparent  green.  The  other  constituents. 


Neptune  Bit 
stable  dye  (among  l 


For  the  duration  of 


lightfastness 


dye  141). 


yet  which  correlates  the  lightfastness  with  its  chemical  composition 


Temperature-Time  History  of  the  Stills 


The  measured  temperature-time  history  of  t 
systems  are  plotted  in  Figures  26-31.  As  can  fc 
figures,  the  general  trend  is  similar  and  thus 


dye-water  system,  the  ambient  temperature  and  the  conduction  of  heat 
from  this  layer  to  the  bottom  ones.  The  glass  temperature  also  follows 


i of  about  130°F  at  1:30  p.m.  (EST)  (Figure  27). 
ire  (Chapter  IV),  this  plate's  purpose  was  to 


greater  than  the  water  surface  temperature  till  about  6:30  p.m.  (EST) 
after  which  it  rapidly  approaches  the  glass  cover  temperature.  Thus 


OJlIJLUOtltUOX 


I 


i 


86 


s plate  thereby 


justifying  the  assumptions  for  equations  (38)  and  (39) . 

The  temperature  of  the  bottom  layer,  in  the  meantime,  continues  to  , 
rise  since  there  is  heat  conducted  from  the  top  layers  during  daytime 
and  also  a small  portion  of  incident  solar  radiation  is  absorbed  by  it. 

both  because  of  rapid  evaporation  of  water  from  the  surface  as  well  as 
the  cooling  of  environment.  At  around  3:30  a.ra.  (EST)  the  bottom  layer 


system  loses  heat  as  a single  mass.  The  intermediate  layer  also 
follows  the  same  trend  as  shown  in  Figure  28. 

The  bottom  layer  temperature  and  the  time  at  which  it  merges  with 
the  top  layer  temperature  is  an  indication  of  the  optical  depth  of  the 

the  Black  dye  (172.3  ppm)  meet  around  1:30  a.m.  (EST)  (Figure  28). 

energy  absorbed  in  different  layers  of  the  dye-water  system.  Figure 
32  shows  the  percentage  of  solar  energy  absorbed  in  different  layers 
at  12:30  p.m.  (EST).  The  Red  dye  (100  ppm)  absorbs  about  75  percent 
of  the  incident  solar  energy  in  the  top  1 inch  layer,  as  compared  to 

same  layer  thickness.  This  results  in  the  top  layer  temperature  in  tin 
Black  dye  system  being  much  higher,  thereby  prolonging  the  time  of 
merging. 


o ooo 


The  above  point  can  be  further  clarified  by  comparing  the 

layer  temperatures  for  control  still  merge  around  8:30  p.m.  whereas  for 
dye-water  system  they  merge  at  around  1:30  a.m.  (EST) . Again,  this 


the  mass  of  the  system  the  top  layer  temperature  drops  less  rapidly  as 
compared  to  that  of  the  dye-water  system.  Figure  33  shows  the 
temperature-depth  profile  for  these  two  systems  for  different  times. 

Figures  34  and  33  show  the  plot  of  temperature-time  hfstory  of 
the  top  layer  temperature  and  glass  cover  temperature  for  dye-water 

cover.  The  higher  top  layer  temperature  is  als< 


u/iwa) 


This  point  will  be  further  clarified  in  the  next 


e temperature  profiles  explained  above  can  be 
te  distillate  output  for  the  two  systems  for  April  15  and  May  15 , 

.s  higher  (about  0.075  lbs/ft2hr)  for  this 
system  as  compared  to  that  from  the  control  still  (about  0.048  Ibs/ft 

he  control  unit  at  night, 


The  distillate  output,  however,  depends  both  on  the  brine  surface 
temperature  and  the  temperature  difference  between  the  brine  surface  and 


3a 

si 


s 

E 


<J>  0 


<t>  Q 


and  q is  given  by  equation  (31) . However,  in  order  Co  calculate  qm  the 

n - hm  A(m.s  - mlg)  (41) 

From  this  plot  the  value  o £ C is  determined  to  be  0.08  for  all  times 


C were  then  subsequently  used  in  the  analytical  model,  whose 
'ill  be  explained  later.  Figure  40  shows  the  experimental  and 
;d  discillate  output  using  the  above  calculated  coefficient  C. 
Once  q^  is  calculated  h0vap  can  be  plotted  against  time  (Figure  38) . 

around  9:30  p.m.  (Figure  39),  a fact  which  is  confirmed  experimentally 
(Figure  37) . 

From  Figure  38  it  is  also  evident  that  hBM  follows  the  trend  of 
surface  layer  temperature  T and  hence  a relationship  between  them  must 


<Jl|-zJj/njaHV  "M3m 


for  various  days  and  different  dyes  is  plotted  against  the  surface 

• »•“*  •»  »•»*  owaVn  «o 

distillate  output  has  been  plotted  against  T for  various  dyes  and  dif- 


artif icially  then  equation  (43)  is  not  valid.  Thus,  tl 


is  changed 


still,  respectively,  for  April  14,  1978.  Each  loss  is  a percentage  of 


o 


)6H9e£n9^9e 


sen  they  are  Identical. 


Such  comparisons  h 


trations  (Figures  46  to  51) . The  data  have  been  plotted  for  completely 
arbitrarily  as  a day  in  which  the  sun  produces  a shadow  for  at  least 


I 


indication  of  how  well  a dye-water  system  absorbs  solar  rj 


. [1  - exp(-a,x*)  ] d A 


radiation  in  the  dye-water  system.  In  Figure  52,  x*  is  arbitrarily 
taken  to  be  1/2  inch.  Thus,  as  x increases  tho  slope  S decreases  and 
the  relationship  is  given  by  the  following  regression  line. 


Figure  53  presents  th 
regression  lines  (Figures 


standard 


1 


found  using  equations  ( 
evident  that  this  incre 


le  distillate  out| 


An  attempt  was  also  made  to  plot  m and  m.  for  different  months 
of  the  year  (see  Figures  54  and  55) . It  is  obvious  that  the  maximum 


increase  are  only  valid  when  the  outputs  from  dye-water  system  are 
compared  from  an  equivalent  control  still  of  depth  10.5  indies. 

In  both  Figures  54  and  55  the  distillate  output  is  for  clear  and  nearly 


- I 


(top-  jj/sqj)  |,u  •oaouijsja 


such  a day  Che  amount  of  radiation  reaching  the  water  surface  is 
is  absorbed  mostly  at  the  surface  of  the  water,  thereby  making  tl 


to  be  expected  since  these  dyes  have  boiling  points  of  about  350#F  and 
above  (413 ; whereas  the  maximum  surface  temperature  recorded  for  dye- 


dyes  to  compare  the  predicted  Cemperature-Cimo  history  of  the  dye-water 
system  with  that  obtained  experimentally.  The  input  to  the  program 


1*1 


by  polynomial  expressions.  These  polynomials  fitted  the  data  well 
C has  already  been  explained.  The  value  of  was  experimentally 

coming  solar  radiation  from  the  water-air  interface.  This  correction 
cent  from  the  water-air  interface  was  obtained.  With  this  analysis 

for  the  whole  day  and  will  obviously  change  during  the  day,  but  for 
simplicity  it  has  been  used  as  such.  The  thermophysical  property 
values  for  water-vapor  viscosity,  density,  etc.  have  been  calculated  at 
an  average  temperature  of  100 6 F (see  Appendix  I for  data  compilation] . 

between  calculated  and  experimental  results,  for  different  dyes,  are 

4®F)  between  the  calculated  and  experimental  temperatures  of  the  top 

simulation  were  average  for  tile  whole  day.  However,  on  certain  days 
there  were  considerable  fluctuations  in  the  wind  speeds  during  day  and 


s _ 
as 


is 


oo 
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There  is  about  4°F  difference  between  the  predicted  and  the  experiments: 
(before  and  after  the  experiments)  was  used. 

surface  temperature  is  slightly  lower  (about  2“F)  than  that  of  layer 
1/2  inch  below  because  of  evaporative  cooling.  There  is  a difference 

pointed  out  that  the  distillate  output  calculated  from  the  computer 


Productivity 


Effect  of  Various  Variables  on  Still 


Once  sufficient  confidence  is  achieved  in  the  validity  of  the 

parameters  like  ambient  temperature,  wind  velocity,  and  various  concen- 
trations of  different  dyes,  on  the  productivity  of  solar  stills. 


The  effect  of  various  ambient  temperatures  (expressed  as  percent 
of  the  actual  ambient  temperature,  T .)  for  March  31,  1978  (Black  dye, 
172.5  ppm)  on  distillate  output  is  plotted  in  Figure  64.  As  can  be 

distillate  output  by  10.2  percent,  while  an  increase  of  Tatiih  by  20 

still.  The  justification  for  this  has  been  provided  by  faulty  argument 

however,  it  should  be  pointed  out  that  the  water  surface  temperature 
does  increase  with  Increase  in  ambient  temperature  but  the  temperature 
difference  butween  the  water  surface  and  the  glass  cover  (the  driving 
potential  for  mass  transfer)  decreases.  Figure  65  shows  Che  plot  of 


surface  temperature  T of  the  water  surface  and  the  difference  between 
the  temperatures  of  surface  and  glass  versus  time  for  two  ambient 
temperatures.  The  difference  in  maximum  surface  temperature  for  these 

EST)  whereas  the  AT(T  - T ) 
greater  productivity)  compared  to  the  higher 


ambient  temperature  conditions.  As  has 

ambient  temperature  condition. 

Effect  of  Wind  Velocity 


The  effect  of  wind  velocity 
cent  with  wind  velocity  increasing  from  0 to  25  mph. 
However,  with  increasing  wind  velocity  (about  25  mph) 

the  effect  of  wind  velocity  on  productivity  [18J. 


e ultimately  reflected  ii 


ce  by  increasing  the  ci 


dye-water  system  on  produc- 
-ll  decreasing  concentration  of 


(top- 
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Che  productivity. 


the  depth  of  the  basin  is  greater  than  15  Inches,  then  the  bottom  layer 


optimum  cost  benefit  ratio  of  using  a particular  dye.  The  problem  is 
one  of  seeking  the  lowest  concentration,  C,  that  will  achieve  a high 
output,  M,  and  it  is  therefore  especially  sensitive  to  the  rate  of 


(51) 


Figure  69  shows  the  plot  of  equation  (51)  with  A equal  to  1.  From  this 
plot  the  concentrations  of  Black,  Green  and  Red  dye  (where  the  value 


concentration  of  Green  dye  over  Black  dye  by  73  percent,  which  is 
directly  reflected  in  the  cost  of  the  dye.  It  should  be  pointed  out, 

valid  for  that  particular  day  (April  14,  1978)  and  will  change  with 
seasons,  with  the  largest  value  of  this  concentration  for  a particular 
dye  being  in  winter  and  lowest  in  summer. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  can  be  di 

1.  By  adding  water-soluble  dyes  tl 
solar  still  can  be  increased,  l 


f increasing 


still  with  no  dye. 

to  be  most  suitable,  both  from  the  point  ol 

Thus,  in  one  month  of  exposure  to  suulight  its  a.  at  0.5  Mm 
exhibited  slight  degradation  and  thus  reduction  in  its  Q,  at 

Based  upon  the  experiments,  a simple  method  of  calculating  the 


(49)  ] - 

An  empirical  relationship  hi 
distillate  output  as  a func 

There  is  no  difference  in  tl 


e distillate  output  from  the 


day. 

An  analytical  model  which  predicts 
excellent  over  the  range  of  parameters  investigated. 

The  effect  of  wind  speeds  on  discillate  output,  as  predicted 


n distillate  oi 


>ut  occurs  with  increasing  con- 


having 


A method  of  predicting  at 
diseiliate  with  concentration  t; 

For  the  spring  conditions  at  Gainesville,  Floric 


i,  respectively. 


It  will  he  worthwhile  to  develop  water  soluble  dyes  which 
have  higher  a.  in  the  wavelength  range  of  0.7  to  0.9  ua 
and  have  higher  X as  compared  to  those  tested  in  the  present 
study.  Besides  higher  x an  attempt  can  also  be  made  to 
develop  dyes  with  greater  lightfastness  so  that  there  is  no 
degradation  by  sunlight. 

A worthwhile  effort  would  be  to  study  the  effects  of  dye- 
water  surfactant  complex  on  evaporation  of  water.  Addition 

with  addition  of  dyes,  there  is  a possibility  of  enhancing 

>uld  also  be  made  to  look  at  different  ways,  both 
physical,  of  extracting  the  dye  from  concentrated 
:o  recycle  it.  Success  in  recycling  the  dyes  may 


Besides  the  use  of  these  dyes  in  solar  distillation  it 
effluents,  frojn  chemical  plants,  by  evaporation, 
liquid  solar  collectors. 


APPENDIX  I 

PROPERTY  VALUES  USED  IN  ANALYTICAL  MODEL 
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APPENDIX  II 

STABILITY  CRITERION  FOR  DIFFERENCE  EQUATIONS 


APPENDIX  II 

STABILITY  CRITERION  TOR  DIFFERENCE  EQUATIONS 
In  Che  calculation  of  the  stability  criterion  for  the  difference 
followed.  The  nodal  equation  for  nth  layer  is: 


Qabs  ° qJitCt)dX  - l",i(t)  exp  ( -a  A Ax*)  dl) 


Subtracting  equation  <11— 1)  from  (34)  yields: 


Rearranging  equation 


tch,  equation  (Ii-5)  results  ii 
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A similar  sore  of  analysis  on  the  glass 


APPENDIX  III 

CALCULATION  OF  ABSORPTION  COEFFICIENT 
DYE  SOLUTION 


APPENDIX 


CALCULATION  OF  ABSORPTION  COEFFICIENT  OF 
DYE  SOLUTION 


TJ 777 


an- 


l-cxp(agt) 

p.  “ spectral  reflectance 
Similarly  the  transmittance  1 __ 


= spectral  absorption  < 


(1  - ag*2  p3  p4) 


(III-*) 


Therefore  the  overall  transmittance  for  the  cell-water  solution  is: 


(III-6) 


<l-o 


(l-o 


thus  knowing  the  values  of  xeM  (empty  cell)  and  lcell/l,yQ 
dye  solution)  the  spectral  absorption  coefficient  a.  can  be 


2,  IFA/Plcnum,  Nt 


Thermophyslcal  Properties  of  Matter.  VoX. 
(1976) . 
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